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The Schwarzschild exterior space-time is conformally related to a direct product space-time, Af 2 x 
S2, where Af2 is a two-dimensional space-time. This direct product structure arises naturally when 
considering the wave equation on the Schwarzschild background. Motivated by this, we establish 
some geometrical results relating to M2 that are useful for black hole physics. We prove that M2 
has the rare property of being a causal domain. Consequently, Synge’s world function and the 
Hadamard form for the Green function on this space-time are well-defined globally. We calculate 
the world function and the van Vleck determinant on Al2 numerically and point out how these 
results will be used to establish global properties of Green functions on the Schwarzschild black hole 
space-time. 


I. INTRODUCTION 


In considering the propagation of scalar waves on the Schwarzschild exterior space-time, two observations have 
significant consequences for any analysis of the problem. 

First is the seemingly mundane observation that one should divide by r. That is, in the usual Schwarzschild 
coordinates, where r is the area radius, technical advantages accrue by considering not the scalar field $ itself, but 
the rescaled quantity dl given by $ = A>/r. (We note that $ satisfies the wave equation 0^$ = 0, where Dg is 
the d’Alembertian operator associted with the metric tensor gap of Schwarzschild.) Following a decomposition into 
spherical harmonic £-modes, the equation for each mode of T* then takes the convenient form of a one-dimensional 
wave equation with an effective potential Vi (and a new radial coordinate r*): 


d^A>i d^^>i 


Vi{r)A!i = 0 . 


( 1 . 1 ) 


Much progress can be made on this basis. See for example Chapter 4 of [T]. 

The second observation is that the exterior Schwarzschild space-time has a non-trivial global causal structure. In 
particular, the exterior Schwarzschild space-time contains circular and other closed null geodesics. This means that 
pairs of spacetime points may be connected by more than one causal geodesic (e.g. a freely falling observer can receive 
a signal emitted in his or her own past). A technical consequence of this is that a large body of the mathematical 
machinery developed for the analysis of waves in curved space-time can be applied only to local (in time) results 
[2]. For example, neither Synge’s world function (a half of the geodesic distance squared) nor the analytic Hadamard 
form [an for the Green function associated with the wave equation are globally defined. Similarly, the existence of 
the photon sphere and trapped null geodesics causes difficulties in the analysis of wave propagation via vector-field 
and energy methods (see Section 4 of W)- 

The main point of the present paper is that considering a geometric interpretation of the first observation (“divide by 
r”) mitigates the difficulties associated in the second observation. The complex causal structure of the Schwarzschild 
exterior is tightly controlled, and we can make significant progress in identifying globally valid versions of the world 
function and the Hadamard form for the Green function. The relevant geometric interpretation is the following: the 
rescaled field 4' is most naturally understood as a solution of a wave equation of the form □§4' + u(a:“)4' = 0 on the 
spacetime with metric g conformally related to the metric of Schwarzschild via 


9ap — 


r ^5a/3, 


( 1 . 2 ) 


and u is a function on the conformal space-time. 

Thus in working with 4*, one is naturally working on the space-time {M,gap)-, where 

M=M2XS2 (1.3) 
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with (in Schwarzschild coordinates) 


■M2 = {(t,r) : t S K,r > 2M}, (1.4) 

M is the mass parameter of the Schwarzschild space-time, S 2 is the unit 2-sphere and the line element associated with 
gap is 


df = r-^ds^ = -^{fdf - r^dr^) + dQ^ (1.5) 

Here, / = 1 — 2 M/r and dfl^ is the standard line element of the unit 2-sphere. We note that this space-time has a 
direct product structure: the metric tensor takes the form 

( 1 - 6 ) 

where (capital Roman indices) and x“ (lower case Roman) are local coordinates on A42 and S 2 respectively. 
Space-time coordinates are obtained by concatenation: x“ = (x^,x“). Likewise, c/ab and 'jab are metric tensors on 
M 2 and S 2 respectively. 

We are thus led to the 2-dimensional space-time that is the focus of this paper. We will refer to this as the 
2 -dimensional conformal Schwarzschild space-time. This is the space-time M2 with line element 

gAsdx^dx^ = dsl = —^{fdt^ — f~^dr^). (1.7) 


Our main result is this: 

Theorem 1 The spacetime {M2,gAB) is a causal domain. 

A causal domain is a geodesically convex domain that is subject to a certain causality condition (see below). A 
geodesically convex domain is a (region of a) space-time where all pairs of points are joined by a unique geodesic. 
Being a causal domain is a technically-advantageous property. For example, Synge’s two-point world function cr(x, x') 
is defined to be one half of the square of the geodesic distance from x to x'. This is well defined only in a geodesically 
convex domain. Similarly, fundamental solutions of wave equations - such as the analytic Hadamard form ma - 
are guaranteed to exist in causal domains (i.e., this is a sufficient, but not necessary, condition for the existence of 
a fundamental solution). Every point x in a space-time has a neighbourhood which is a geodesically convex domain 
{normal neighbourhood), and so a is always defined locally. For renormalization within quantum field theory on 
curved space-time, coincidence limits of two-point functions play a central role: thus local existence is sufficient [7]. 
However, in the study of wave propagation, global existence of a and consequently of the Hadamard form is extremely 
desirable. Such existence holds in a causal domain. 

The property of a space-time of being a causal domain seems to be an uncommon one. In particular, Schwarzschild 
space-time is not a causal domain, and so the world function and the Hadamard form are only well-defined in normal 
neighbourhoods. As a consequence, one must resort to other techniques in order to calculate the Green function 
outside a normal neighbourhood in Schwarzschild (e.g., mi[)- 

The main result of this paper is the proof of Theorem 1. In addition, we study the geodesics of AI 2 and numerically 
calculate the world-function and the van Vleck determinant, which are well-defined globally on Al 2 - The van Vleck 
determinant is an important bitensor appearing in the Hadamard form for the Green function which gives a measure 
of the focusing of the spray of null geodesics emitted from a point m- For the Green function presentation and van 
Vleck determinant results, we focus on the case of a massless scalar field. In a separate paper we use these results in 
order to obtain a globally-valid expression for the Green function in Schwarzschild space-time and derive some of its 
properties m- 

In Sec|H we give details of how the space-time M 2 arises naturally when we consider the wave equation of the 
exterior Schwarzschild space-time, and we describe associated Green functions. In Sec |HI| we recall the definition of 
a causal domain, we determine properties of the geodesics on AI 2 and we use these to give a proof of Theorem 1. 
In Sec |IV| we provide the numerical calculation of a and the van Vleck determinant in M 2 . We then make some 
concluding comments. 

Throughout this paper we choose geometric units c = G = 1 and metric signature (—h -I--I-). 
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II. FROM THE WAVE EQUATION ON SCHWARZSCHILD TO M2 


Our main aim in this section is to show how the space-time M 2 arises naturally in the study of the scalar wave 
equation on the exterior Schwarzschild space-time. The line-element for this space-time can be written as 


ds^ = — dri) r'^dfll, 


( 2 . 1 ) 


where, as above, M is the mass of the black hole, / = 1 — 2M/r and is the standard line element of the unit 

2-sphere. The tortoise coordinate r* is related to the area radius r by 


dr 

dr^ 


= /• 


( 2 . 2 ) 


A massless scalar field d) propagating on Schwarzschild space-time satisfies the Klein-Gordon equation 

= 0 . 


(2.3) 


Our principal interest is the study of Green functions for this equation, continuing a line of research that approaches 
the self-force problem through the study of Green functions on black hole space-times (e.g., [HIEHlEllE] ) (although 
the global behavior of Green functions also finds applications in other areas, such as in the calculation of the response 
of a quantum “particle detector” Cl)- For this reason, we will immediately focus on Green functions rather than on 
the scalar field itself. 

The retarded Green function Gii{x,x') is a solution of the inhomogeneous version of ( |2.3| ) with a Dirac-delta 
distribution as a source, satisfying certain causal boundary conditions: Gr(x,x') vanishes ouside the past light cone 
of X. In the coordinates above, the wave equation for Gr reads 


^ j-iQ^Gr 


dt^ 


9r* 


2 OGr 1 ^ 2 ^ ^TT A Al 

" -f -^\/^Gr = — ^S2ix^ - x^ ) 


dr A, 


j.2f 


S 2 (x“,x“ ), 


(2.4) 


where is the Laplacian operator on the unit 2-sphere, x^ = {t,ri,) are coordinates on the Lorentzian 2-space 
(i.e. the 2-space that arises by factoring the 4-dimensional space-time by the action of the SO{3) that generates the 
spherical symmetry) and x“ = (0, 4>) are coordinates on the unit 2-sphere. 

It is standard at this point to rescale the scalar field $ by a factor r: this removes the first order derivative from the 
wave equation (2.4) and simplifies the use of a variety of methods from the theory of Partial Differential Equations 
(e.g. WKB approximation methods, frequency domain analysis)^. It is worth considering the geometrical significance 
of this step. The appropriate rescaling amounts to making a conformal transformation of the metric and considering 
4' to be the solution of a corresponding equation on the conformal space-time. That is. 


□,4> = 0 


□e4'-5^^ = o, 


(2.5) 


where gajs = r '^ga/ 3 - Thus we may consider that the evolution takes place on the spacetime with line element (1.5). 
By general properties of Green functions in conformally related space-times, we can write laE! 


Gr= - jGr{x,x'), 

r . IT'' 


( 2 . 6 ) 


where Gr{x, x') is the retarded Green function for the conformally invariant wave equation on the space-time with line 
element (Note that since the Ricci scalar R vanishes in Schwarzschild space-time, the ordinary wave operator 

□g and the conformally invariant wave operator Dg — R/6 coincide. Moving to the conformal space-time naturally 
gives rise to the conformally invariant wave equation on that space-time - which is the second equation of ( |2.5| ).) We 
will refer to this as the conformal Schwarzschild space-time. Using this rescaling, we find 

“ ^7)^ + 77 % + “ —)Gr = -47r(52(x^ - x'^')(5s2(x“,x“'). (2.7) 

ot^ r 


^ Note that we could rescale $(x) by any constant multiple of r. The corresponding rescaling of Gji{x,x') likewise involves the factor r, 
multipl ied b y any term that is constant with respect to the point x. For reasons of symmetry, we choose this ‘constant’ multiple to be 
r': see (|2.6||. 
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In addition to rendering the wave equation more tractable, the conformal rescaling introduces a very useful sim¬ 
plification of Synge’s world function. The world function between two space-time points is half the square of the 
geodesic distance along a specific geodesic joining the two points mg. As pointed out earlier, this is well-defined 
only on regions of space-time in whi ch t here exists a unique geodesic between each pair of points. The direct product 
structure of the metric induced by (1.51 yields 


CT4 = a{x^,x^') + ^ 7 ^ 


( 2 . 8 ) 


where is the world function of the conformal Schwarzschild space-time and a{x^,x^ ) is the world-function of 


the 2-dimensional Lorentzian space-time with line element (1.7). This is the 2-dimensional conformal Schwarzschild 
spacetime, denoted by M. 2 - We note also that in (2.8), 7 is the geodesic distance on the unit 2-sphere. Then, by 
conformal invariance of null geodesics, a null geodesic connects x, x' in Schwarzschild space-time if and only if a null 
geodesic connects the corresponding points of the conformal Schwarzschild space-time. 

Next, we examine the space-time geometry of M. 2 - We derive the geodesic equations and some of their properties 
and use these to prove that AI 2 is a causal domain. 


III. GEODESICS ON M2 AND THE PROOF OF THEOREM 1 

This section gives the main results of the paper. We begin by recalling the definition of a causal domain and 
showing that M .2 satisfies the first of the two relevant criteria. We then write down the geodesic equations of A 42 
and describe briefly, in terms of these equations, what it is we need to prove in order to establish that M .2 is indeed 
a causal domain. We then give the relevant proof, establishing the required properties of spacelike, null and timelike 
geodesics in four separate subsections. We note that the argument in the case of timelike geodesics is considerably 
longer than that for spacelike and null geodesics: issues of uniqueness and existence are dealt with separately. 


A. Causal domains 


As noted in the Secim the world function of the conformal Schwarzschild space-time decomposes according to 


the 2-dimensional conformal space M 2 - We will argue here that M 2 a causal domain [2], so that a is defined 
globally. This yields the major advantage that the Hadamard form for Green functions on M 2 (see Eq.( |4.5[ ) below) 
is globally valid: this in turn will yield a globally valid expression for the retarded Green function on Schwarzschild 
space-time [13]. 

A spacetime {M^gap) is a causal domain if 
(GD-i) any two points p,q of M are joined by a unique geodesic and 






2.8). This form is of course only locally valid in the 4-d space-time, and involves the world function a{x^,x'^ ) of 


(GD-ii) for all pairs of points p,q in M, the set J^{q) H J {p) is a compact subset of Al, or is empty. 

The set J~^{q) is the future emission of q G M, defined to be the closure of the chronological future D^{q) of q: 


J+{q)=D+{q), 


= {p & M : there exists a future-directed time-like geodesic from q to p}. 


(3.1) 


The past emission of q £ M, J~{q), is defined analogously to T+(g), but with ‘future-directed’ 
‘past-directed’. 

The space-time of interest is AI 2 = (I^^jffAs) where in coordinates x^ = (t, r*). 


in (3.1) replaced by 


9ab = 




(3.2) 


Then M 2 is globally conformally flat: gAB = f ‘^fVAB for all space-time points, where tjab is the metric of 2- 
dimensional Minkowski space-time. It is then straightforward to show that 

= {(i>c) G : k* - c(g)| < t - t{q),t> t{q)}, (3.3) 


and consequently, for any pair of points p,q G M 2 , J'^{q) H J~(jp) is either a closed diamond, a single point, or is 
empty. Thus this set is always either compact or empty, and so the condition (GD-ii) for M 2 to be a causal domain 
is satisfied. 

It is less straightforward to show that any two points of M 2 are connected by a unique geodesic, but we will prove 
that this is the case. 
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B. The geodesic equations 


In order to study the geodesics on AI 2 , we rewrite the line element as 

dsl =-A{r*){dt^ - drl), A{r^) = 


(3.4) 


In (2.2), we choose a constant of integration so that the tortoise coordinate satisfies r* = 0 at r = 3M. Thus 

/r- 2 M' 

r» = r — 3M + 2M In 


M 


(3.5) 


and we note that r* —>■ —00 as r —>■ 2M+, r* —>■ +00 as r —>• + 00 . Then > 0 for all r, S M, is increasing on (— 00 ,0), 
is decreasing on (0, + 00 ) and has a global maximum of I/27M^ at r» = 0 (r = 3M) - see the plot in Fig{^ The radii 
r = 2M and 3M, of course, respectively correspond to the radii of the event horizon and the unstable photon orbit 
in Schwarzschild space-time. 

The geodesic equations are (with an overdot for derivative with respect to an affine parameter s along the geodesic, 
and a prime for derivative with respect to argument) 


•• . 2 

r* -b - 


with the first integral 


o' 

II 

(3.6) 

A' 

(3.7) 


rj) = e. 

(3.8) 


and where e = 0, —I, -bl for null, timelike and spacelike geodesics respectively. 

Our aim is to prove that give n an y pai r of points (fi,r»i), (i 2 ,?'* 2 ) G there exists a unique value of e and a 
unique solution of the equations (3.6), (3.7) that connects these two points. Without loss of generality, we take ti < ^ 2 - 


The equation (3.6) may be integrated to yield 


At = E, 


(3.9) 


where E is constant along the corresponding geodesic. (This conserved quantity corresponds to the existence of the 
Killing vector ^ in M. 2 -) Since A > 0, this means that t is monotone along any geodesic. Without loss of generality, 
we choose E > 0,so that t > 0 for all geodesics. This amounts to choosing s and t to be co-synchronous along geodesics 
with E ^ 0. Furthermore, since E is constant, either i = 0 everywhere along the geodesic, or t > 0 everywhere. 

In the case that t > 0 everywhere along the geodesics in question (we consider the special case t = 0 in the next 
subsection), we can reduce the problem to one dimension, as we can take t to be the parameter along the geodesic 
and reduce the geodesic equations to the single equation that arises from (3.8): 


dr* \ A 


(3.10) 


The remaining geodesic equation may be taken to be t = E/A: given a solution of (3.10), this simply serves to 
determine the parametrisation t = t{s). 

In the general case (i.e. t > 0), two points q = (fi,r*i) and p = (f 2 T* 2 ) in A 42 are null separated if \t 2 — ti| = 
|r *2 — r*iI, timelike separated if 1^2 — | > |r *2 ~ ^*11 and spacelike separated if |t 2 — | < |r* y — r* i |: notice that any 

given pair q,p has a unique separation character^. For geodesics with t > 0, it is clear from (3.10) that 


dr* 

dt 



e 



(3.11) 


^ In a general space-time, a pair of points p,q are null (timelike, spacelike) separated if there is an everywhere null (timelike, spacelike) 
curve that connects the points. Pairs of points do not necessarily have a unique separation character: consider e.g. points (t, r) = (ti, 3M) 
and (f,r) = (t2,3M) on the photon sphere r = 3M in Schwarzschild space-time, with fi,t 2 chosen so that a circular null geodesic at 
r = 3M passes through both points. These are also connected by a time like curve with r = 3M and 0,(f) constant. However, since 
24(r*) > 0 for all r* E M, it is straightforward to show that the separation character is unique in j\42 and that it may be described as 
stated in the text. 














6 


It is also clear from Eq.(3.8) that the case t = 0 and r* 7 ^ 0 corresponds to e > 0, i.e., a spacelike geodesic. 


It follows that null (respectively timelike; spacelike) separated points can be connected only by null (respectively 
timelike; spacelike) geodesics. This elementary observation allow us to decompose the proof that (CD-i) holds in a 
useful way, by considering the pair of points pi = (ti,r*i),p 2 = (i 2 )?’* 2 ) from the point of view of their separation 
character. That is we prove that if \t 2 — ti| = \r ^2 ~ ’’♦il) then there is a unique null geodesic from pi to p 2 ', if 
\t 2 — ti| < |r »2 — r*i|, then there is a unique spacelike geodesic from pi to p 2 and finally if \t 2 — ti| > |r *2 — then 
there is a unique timelike geodesic from pi to P 2 - This exhausts all possibilities and so establishes (CD-i). 

Before proceeding, we note that we can relate the radial coordinate r and the affine parameter s by combining 
Eqs.(|3.9) and (|3.8|) to obtain 


-^=E^ + eA, 

with A = A(r*) acting as a radial potential. We plot this potential in Fig[^ 


(3.12) 





r/M 


FIG. 1: Radial potential A in Eq.(3.12| as a fnnction of r/M. Its maximum is at r = 3M. 


C. Null separations - existence and uniqueness 


It is trivial to show that null separated pairs of points are connected by a unique null geodesic. The required 
geodesic is 


I ^*2 ^*1 /.A J. \ 

rAt) = r*i -I- —- 7 — (t - ti). 

ty — ti 


(3.13) 


D. Spacelike separations - existence and uniqueness 


We have two cases to deal with here: ^2 = ti and t 2 > ti. In the former case, we must have i = 0 everywhere along 
the geodesic (since t > 0). Then e = -1-1 in (3.8 1 and we can integrate to obtain the general form 

r(s) = M(1-I-cosh(s — So))i (3-14) 

where s, the affine parameter along the geodesic, represents (non-dimensionalized) proper distance. Initial and 
terminal values of s then yield the unique geodesic connecting (ti,ri) and (ti,r 2 ). 

For the case f > 0, we may take r '*2 > without loss of generality. Then we seek a unique value of E and a 
unique solution of 


dr* / A 


(3.15) 


that connects these points. For a given value of E and a given initial point q = (ti,r*i), we can formally write down 
the solution of (3.15): 


t2{E) — ti — 


rr.5 
Jr, I 


E 


AE^ + A{rA 


dr,.. 


(3.16) 
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As emphasised, this equation yields the value of t 2 as a function of E with and r *2 taken to be fixed. We can 

then calculate 


dt2 

dE 



A(r*) 

(i;2 + A(r*))3/2 


dr^ > 0 , 


(3.17) 


so that the function E i— t 2 {E) is monotone increasing. We note also that ^ 2 ( 0 ) = ti. For any given pair of values 
r*i,r* 2 , the function A(r*) has strictly positive lower and upper bounds on [r*i,r* 2 ]. It follows that 


lim ^ 

E ^+00 + A(r,) 


= 1 


(3.18) 


uniformly on [r»i,r* 2 ], and so for r*i,r *2 fixed, 


lim t 2 (E) — ti 

E^+00 


|•r,2 

lim / 

E—>.+oo 





dr* 


r*2 -r*i. 


(3.19) 


It follows that E H. t 2 {E) — ti is a one-to-one mapping of (0, 00 ) onto ( 0 ,r *2 — r*i). Therefore given any pair of 
point s satisfying the spacelike separation condition t 2 — ti < r *2 — r*i, there exists a unique value of E such that 
(3.16) provides the unique spacelike geodesic connecting these two points. 


E. Timelike separations - general issues and uniqueness 


Let us now consider the case of timelike geodesics. This case is more complex than the others, as the right hand 
side of (3.10) can vanish along the geodesic, depending on the value of E. Recalling that A{r*) has a global maximum 
of 


1 

27M2 


at r* = 0, we define the critical value of E to be 


Er = 


3V5M' 


(3.20) 


This gives rise to the following classification of timelike geodesics: 

(i) Supercritical timelike geodesics. These have E > Ec, and so - from (bdol - > 0 for all t. Hence 

t !->• r*{t) is monotone along these geodesics. As the right-hand side of ( |3.10p is uniformly bounded and is 
bounded away from zero, these geodesics are defined for all t G M and satisfy either lim(_,.±oo r* (t) = ±00 or 
limt_,.±oo r»(t) = Too depending on the sign of The supercritical geodesics overcome the potential barrier 
at r* = 0 . 


(ii) Critical timelike geodesics. These have E = E^, and (3.10|) may be written in the form 


dr*'\^ _ _ A{r*) 
dt ) A( 0 ) 


(3.21) 


The unique geodesic for which r*{to) = 0 at some to G ® is the static geodesic r*{t) = 0. Thus all other critical 
geodesics are confined to either < 0 or r* > 0. These satisfy either limt_>_oo r*{t) = 0 or limt_,.+oo r*{t) = 0. 


(hi) Subcritical timelike geodesics. These have E < Ec- For a given E G (0,i?c), there exist unique values 
r*^{E), with r*_|_(if) > 0 and r*_{E) < 0, f or which 1 — A{r*±)/E‘^ = 0. These values correspond to turning 
points of the geodesics: differentiating (3.10) and noting that r'^{t) ^ 0 yields 


d'^r* _ A'{r*) 

dE ~ 2E^ 


(3.22) 


so that subcritical geodesics in r, > 0 have a global minimum at r*_|_(£’) > 0 , and subcritical geodesics in r* < 0 

have a global maximum at r*_{E) < 0. These geodesics reflect off the potential barrier at r*^. We note also 

that for a subcritical geodesic with r* > 0 (respectively, r* < 0 ), there is a unique value to such that r'*{t) < 0 

for t < to, r*{to) = r*_f_ and r(.(t) > 0 for t > to, and r*{t) —>■ -l-oo as t —)■ ±00 (respectively, r(.(t) > 0 for t < to, 

r*{to) = r*_ and r*(t) < 0 for t > to, and r*{t) —>■ —00 as t —)■ ± 00 ). 




















The simple observation that 


0 < El < E 2 


Mr*) _ Mr*) 
El El 


(3.23) 


allows us to rule out multiple crossings, and thus demonstrate uniqueness, for many cases of pairs of geodesics. For 
example, for a pair of supercritical geodesics with E^ < Ei < E 2 , both satisfying ^ > 0 , this shows that (with 


obvious notation) 


dr^ 

dt 


< 


dr^ 

dt 


- at a point of intersection. Thus r* = r*]^(t) always crosses r* = r* 2 (t) from above, and 


so the geodesics meet at most once. Variations on this argument immediately rule out multiple intersections except 
for the case of a pair of subcritical geodesics. 

We will not go through all these possible variations, but by way of example, deal with one case, that of a super¬ 
critical or critical timelike geodesics meeting a subcritical timelike geodesic twice. For example, consider a critical or 
supercritical geodesi c (so that E = E 2 > Ec) with ^ < 0 and a subcritical geodesic confined to r* > 0 (and with 
E = El < Ec). By (3.23), the supercritical geodesics either (a) crosses the subcritical geodesic from above at some 
r* < r* + , the minimum of r* on the subcritical geodesic: a second crossing could only occur on the decreasing branch 
of the subcritical geodesic, but this is ruled out by comparison of slopes; or (b) first crosses the subcritical geodesic 
on its increasing branch at some r* > r*_|_, and no other crossing is possible. 

Thus, to establish uniqueness, it remains to consider the case of two subcritical geodesics intersecting twice. This 
corresponds to a pair of particles Pa and Pb, moving on timelike geodesics 7 ^ and 73 with energies Ea < Eb < Ec, 
which meet at the points (r*]^,ti) and (r* 2 ,t 2 ) with ^2 > fy- We argue as follows that this is not possible. 

The argument is quite straightforward, but has some technicalities (relating to the second and third points listed 
below) that are detailed below. We consider the case of two subcritical timelike geodesics confined by the potential 
barrier to r* > 0. The case r* < 0 is similar. 

The uniqueness argument is based on the following observations: 


0-1. The relation between the energies and slopes of the trajectories at a point of intersection is constrained by 

0-2. The motion of the particle is symmetric about the minimum. That is, r*(t+ — t) = r^{t+ + t) for all t G M, 

where the minimum r* = occurs at t = fy. We refer to as the arrival time. 

0-3. The minimum value of r, along a subcritical geodesic is a decreasing function of E. That is, rJ_^_{E) < 0. See 

- recalling that ^ > 0 . 

0-4. Consider all subcritical geodesics, indexed by P < Ec, which decrease through a point (ti, r*i). Define Ai{E) = 
t+M) ~ fyi where t^{E) is the arrival time of geodesic E, i.e. the time at which geodesic E meets its minimum. 
Then Aj(i?) > 0, so that lower energy particles meet their minimum before higher energy particles emitted 
from the same point at the same time. Likewise, subcritical geodesics which increase through a point (t 2 ,?’* 2 ) 
satisfy A^E) > 0 , where A 2 M) = fy — t+{M- 

Subject to these conditions, there are then four possible ways in which the two subcritical geodesics can meet twice, 
illustrated in Figure 

We rule out these four cases as follows. 



(a) The geodesics intersect at (ti,r*i) and (t 2 ,?'* 2 ) with < 0 < t 2 - By 0-4 above, 

Ai{Ea) = ti_{EA) — ti < t+^EB.) —ti= Ai{Eb), (3.24) 

since Ea < Eb (this is correctly represented in the diagram). But also 

A2{Ea) = fy — t+{EA) < fy — t+{EB) = A2 {Eb). (3.25) 

Adding the inequalities leads to the contradiction fy — ti < ^2 — ti, and so case (a) is ruled out. 

(b) Case (b) also contradicts 0-4, as 75 - the higher energy geodesic - reaches its minimum before ^a does, having 
started from a common point. 

(c) In case (c), the second intersection point violates the slope inequality of 0-1: on an increasing branch, the higher 
energy geodesic 75 can only meet the lower energy geodesic ^a from below. 

(d) The same argument as case (c) applies to rule out case (d). 
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FIG. 2 : The figure shows the four inequivalent ways in which a pair of subcritical timelike geodesics could intersect. The 
diagrams are representative - our aim is to prove that such intersections cannot occur - but note that 0-2 is respected in the 
diagrams. In each graph, r* runs along the vertical axis and t on the horizontal, and Ea < Eb in each case {Ei being the value 
of E along the geodesic 7^. Notice that this gives rise to a lower minimum for 'ys- By a translation, we fix the minimum of 'ja 
to be at t = 0. 


To complete the argument, it remains to establish the claims of 0-2, 0-3 and 0-4. 

The first of these is trivial: each branch of the geodesic reaching/emerging from the minimum is the unique solution 
of the relevant time-reversed initial value problem. The condition (0-3) follows immediately from (3.30) below: 
establishing (0-4) requires more effort. 

So consider a subcritical geodesic with energy E < Ec (recall that Ec is the maximum value of the energy E 
such that the timelike geodesic bounces off the potential) with initial value = r*(ti) > 0, and which is initially 
decreasing, so that 


dr* 

dt 


= -\ 1 - 


A(r*i) 


£2 


(3.26) 


That is, the subcritical geodesic starts off at an initial point (ti,ri) and bounces off the potential barrier at the 
periapsis r* = r*_|_. For this case we must have -y/A{r^i) < E < Ec- 

The value of r* along the geodesic decreases monotonically until it reaches its global minimum r* = r*_|_(if) at time 
t = t+^E). A formal integration yields 


Edr^ 


or, returning for convenience to the area radius coordinate r 

rri 


Ai{E) = t+{E)-ti= f ' ,_, 

Jr,+{E) \/A2 - A(r*) 

I 

^ 1 dr 


Here, r_|_(if) is the unique root of 


-r + 2M = 0 


(3.27) 


(3.28) 


(3.29) 


with r > 3M. With E < Ec, (3.29) has three real roots; one negative root, denoted r„, and two positive roots: 
r_ G (2M, 3M) a nd r^ G (3M, -Poop We note that the area radii r± correspond to the tortoise radii r*±. Implicit 
differentiation of (3.29) yields 


= -Erl 

dE + 


r+ J 


-1 


< 0 . 


(3.30) 
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This inequality has a simple physical interpretation: subcritical particles with higher energy penetrate deeper into 
the potential barrier. Note that this establishes 0-3. 


Monotonicity allows us to consider t_|_ as a function of r_|_, and substituting for from (3.29) yields 

l-ri ^ 3 / 2^2 




dr 




ir-r+){r^+r+r-^^) 


1 1/2 ’ 


(3.31) 


where <5i(r_|_) = Ai(E). As the particle increases away from the minimum, it reaches r *2 at time t 2 > > ti and we 

can write 


A2{E) — t2 — t+{E) — 
or equivalently (with i52(r_|.) = A 2 {E)) 

fr2 

S2{r+) = I 


Edr^ 


r.+ - A(r*)’ 


rV^l-my^dr 


J( \( 2 \ 2Mr^^ \ ^ 

+ [(r - r+)(r^ r+r - ) 

To establish 0-4 above, it remains to show that for any fixed xq> x > 3M, the function 


r 

X i-> St{x) = / 

J X 


^3/2{i-m)-yr 


{T-x){r^+xr-^) 


1/2 


is decreasing. This follows from the fact that (taking Xg = ri and x = r+) 


and that < 0, as we know from (3.30). The case for A 2 {E) is identical. 

By rescaling, and without loss of generality, we set M = 1 for the remainder of this subsection. If we introduce 


(3.32) 


(3.33) 


(3.34) 


(3.35) 


X = r+, 

then considering the nature of the roots of the cubic equation ( 3.29| ) allows us to write 

r_ = b{x)^ Tn = —X — b(x). 

where 


bix) = -l-l + 


a; -f 6 
X — 2 


We note that 


Then we can write 


St{x) = 


b{x) < X for all a; > 3. 

r^/'^dr 


(r — 2)i/r — Xy/r — b{x)y/r -|- x -|- b{x) 


(3.36) 

(3.37) 

(3.38) 

(3.39) 

(3.40) 


Thus to prove that A42 is a causal domain, we must show that (5 is a decreasing function of x on the interval (3, a;o), 
for all a;o > 3. 

Let us recall the meaning of 5t{x). For each fixed xg, this measures the elapse of coordinate time t as a, particle 
falls from r = a;o to its minimum r = a; = r_|_ (which by symmetry, is the same as the time to reflect off the potential 
barrier at r = a: and move outwards to r = Xg). Instead of dealing with the coordinate time t, we shall find it easier 
to deal with the proper time r. 
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It is straightforward to show that the function corresponding to 6t is 


6r{x) = 


1 


dr 


^{x) _l_ a; _l_ 


(3.41) 


and the problem now becomes to prove that x i—>■ Sr{x) is a decreasing function of x on (3,a;o) for all Xq > 3. Note 
that we have emphasized that E' is a function of x. Indeed since r = x = r+ solves (3.29), we have 


£;(x) = 


x-2My/^ 

x^ ) 


(3.42) 


The advantage of moving to the proper time representation is that Sr may be written in terms of a familiar special 
function (see Eq. (3.147-8) of [T5] : 


Sr{x) = aF[ip,q) 


where 


= 1 7r 


dO 


qsm 


is Legendre’s incomplete elliptic integral of the first kind m ^ and 


a = 


_ /2b(x + b)\ 


1/2 


Ey/ir- - r„)r+ \ x + 2b J 

Ip = sin“^ (p> 0, (p= 

^ (r_ - r„)(xo - r+) _ {x -f 26)(xp - x) 

^ (r+- r„)(xo - r_) (2x-b &)(xo - &) ’ 

^ r_(r+ - r„) _ b[2x F b) 
r'+(r_—r„) x(x-b 26) 

It is straightforward to establish that 

0<(/J<l, 0 < Ip < 7r/2, 0 < q < 1 for all Xo > X > 3. 


(3.43) 

(3.44) 

(3.45) 

(3.46) 

(3.47) 

(3.48) 

(3.49) 


We shall now use the exact result in Eq.(3.43) to prove that Sr is a decreasing function of its argument, and so 
prove 0-4. 


We note the partial derivatives 

dF{ip, q) 
dtp 

dF{ip, q) 
dq 


where 


1 


\/l — q sin^ Tp 
1 


2q{l - q) 


E{ip,q) - (1 - q)F{ip,q) - q 


cos Ip sin Ip 
\/l — q sin^ Ip 


E{tp, q) = [ \Jl — q sin^ 9d9 

Jo 


(3.50) 

(3.51) 

(3.52) 


is Legendre’s incomplete elliptic integral of the first kind. Then we find, using (3.46), 

a 


S'rix) = 


2(1 - g)Vl - f 


((1 - q)(p' - q'(p{l - (fi)) 


+ {a - 


a q 

2 i 


F{'ip,q) 


aq 


2q{l - q) 


E{ip,q). 


(3.53) 


^ Here and below, we follow MATHEMATICA’s |16| definition of the elliptic integrals, which differs slightly from that in |15| . 
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To establish the decreasing nature of Sr, we show that 


Ti = (1 - q)(p' - q'ip(l - (p) 


and 


T 2 = PF{tp,q) +-fE{tlj,q) 

are both negative for all x G (3,xo) and all xq > 3. Note that we have introduced 

f / 

aq 


P = a - 

2 q 


7 = 


2g(l - q)' 


For ease of reading, we separate these results into two lemmas. 
Lemma 1 7i < 0 for all x G (3,xo) and all Xq > 3. 

Proof: A straightforward calculation shows that 

— 16x0 


Ti = 


(x - 2)3 Jf^(3 + J^Y{2h - 2xo)2 


X (u + 


X + 6 

X — 2 


v), 


(3.54) 

(3.55) 


(3.56) 


(3.57) 


where 


(3.58) 

(3.59) 


u = {x — 3)^ + 5(x — 3)^ + 5(xo — 3)(x — 3) + 9(xo — 3), 

V = (x — 2)((x — 3)^ — 3(xo — 3)). 

It is immediate that u > 0 for all x G (3, xq) and all xq > 3. If (x — 3)^ > 3(xo — 3), then v is non-negative, and 7i is 
negative, proving the lemma. If (x — 3)^ < 3(xo — 3), then 7i is negative if and only if 


X + 6 2 / 2 : -I- 6 \ 2 

u>-\ - -V u >1 - -\v. 


X — 2 V X — 2 . 

Expanding terms shows that this last inequality is equivalent to 

(34 -h 16(xo - 3))(x - S)'^ -h 128(xo - 3)(x - 3)^ -h (16(xo - 3)^ -f 144(xo - 3))(x - 3)^ > 0, 
which is clearly true for all x G (3,xo) and all xq > 3. This completes the proof. 

Lemma 2 T 2 < 0 for all x G (3, xq) and all xq > 3. 


Proof: From the dehnitions (3.44) and (3.52), we can write 

(3 + "/ — ygsin^ 9 


T2 = 




de. 


gsm 


We hnd that 


and so 


79 = 


aq 


2(1 - q) 


< 0 for all X > 3, 


' + 7 ~ 79sin^ 9 < /3 + j — yqsin^ if, 9 G (0, ip). 


Using (3.46), this allows us to write 


72 < (/3 + 7(1 - qp))F{ip, q), 

and so the lemma is proven if we can establish the inequality 


/3 + j{l-qip) = - 




1-9 


< 0 for all X G (3, xq), xq > 3. 


(3.60) 


(3.61) 


(3.62) 

(3.63) 

(3.64) 

(3.65) 

(3.66) 
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We find that 


where 


(a7 

^2 


Vl“9/ {x-2YJ^h{x + b){2x + b){xo-h) M x-2 


y = (xo + l)(a; - 3)^ + 5xo(a; - 3) + 3(4xo - 3), 
z = —((xo — 3)(x — 3) — 9)(x — 2). 

We see immediately that y > 0 for all x € (3, Xq), Xq > 3. Thus 72 < 0 if z > 0. If z < 0, then 


y + 



z > 0 


<S4> 


> 


/x + 6 
\x-2 


(3.67) 


(3.68) 

(3.69) 


(3.70) 


Expanding terms shows that this last inequality is equivalent to 

(xo-l)(x-3)^ + (llxo-18)(x-3)3 + ((xo-3)2 + 12(xo-3)+9)(x-3)2+3(xo-3)(5xo + 18)(x-3)+9(xo-3)(2xo-3) > 0 

(3.71) 

which is clearly true for all x € (3,xo) and all xq > 3. This completes the proof of the lemma. I 

This completes the proof that timelike separated pairs of points in AI 2 are connected by at most one timelike 
geodesic. We complete the proof of Theorem 1 in the following subsection by proving that timelike separated pairs of 
points are connected by at least one timelike geodesic. 


F. Timelike separations - existence 


We consider two points (ti,r*i), {t 21 X 1 , 2 ) which satisfy the time-like separation condition 

|r*2 - r„i\ <\t2-ti\=t2- h, 


(3.72) 


where we impose, without loss of generality, ti < t 2 - Our aim is to prove that there exists a value E-^ ^ oi E and a 
solution of 


yl(r*) 

dt ) ~ E 2 


(3.73) 


which satisfies r*(ti) = r*i and r*(f 2 ) = r* 2 . In addition, we specify that r*i < 0 (corresponding to r{ti) < 3M): the 
arguments carry over, mutatis mutandi, to the case r*i > 0 . 

We introduce r* = r*(t), the solution of the initial value problem (IVP) 


/_ _ A(r«) 

\ dt ) 


r*(fi) = r*i. 


(3.74) 


That is, this geodesic has energy E„ = yjA{r„i). This geodesic exists and is uniquely defined on K. It has the 
distinguishing property that r* = f„{t) has a global maximum of r* = r*i which is attained at t = ti. Thus this 
geodesic is monotone decreasing on (ti,-|-oo) and so r*(t 2 ) < J"*!- We note that r* is also the unique solution on 
[ti, -foo) of the IVP 


dt'^ 
r* (ti) 


A'{r„) 

E 2 


' *ii 


dr„ 

dt 


= 0 . 


(3.75) 


The initial condition r(,(ti) = 0 fixes E"^ = ^(r*i). 

We establish existence by considering different values of r *2 as delimited by f*(t 2 ) and by r*i ± (t 2 — ^i)- 
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(i) r^{t 2 ) < r *2 < + {t 2 — ti). The key here is that the geodesic we seek is initially non-decreasing. Consider 

the IVP 


dPr^ 

dt^ 

r* (h) 


A'{r,) 

dr^ 



(3.76) 


For each E G [iil*,-|-oo), a solution of this IVP is uniquely defined on [ti,-|-oo). Further more, this solution 
depends continuously on the parameter E, subject to dependence of the coefficients of (3.761 on E, which 
holds for E G [if*, -l-oo) (see e.g. Section 2.3 of [H]). That is, for t 2 > ti, the mapping 


E : [if*, -l-oo) (r*i - (^2 - ii),?'*i -f (t 2 - h)) : E i-G r*(t 2 ; E) 


is continuous, where r*(t 2 ;if) is the solution of (3.76) evaluated at i = ^ 2 - Th e cod omain of this mapping is 
determined by the null geodesic envelope. In the limit E —>• -foo, the solution of ( 3.76[ ) approaches the outgoing 
null geodesic through (ii, r*i), and we have limB_>+oo F*(i 2 ; E) = r*i -I- (^2 — ii)- Similarly, we have r~*(i 2 ; if*) = 
r*(i 2 )- By the continuous dependence result mentioned above, it follows that for all r *2 G (r*(t 2 ), ^*1 -I- (^2 —ii)), 
there exists E 2 G (if*,-|-oo) such that r*(< 2 ;if 2 ) = ^* 2 - This gives use the required timelike geodesic connecting 
(ti,r*i) and (t 2 ,r* 2 )- 


(ii) r*i — (^2 — ii) < r *2 < r*(i 2 )- Existence is proven in this case by applying the argument of case (i) to the IVP 


d^r* 

dE 

r*(ti) 


A'(r*) 

E^ 

dr^ 



(3.77) 


So we conclude the following: Given any pair of events (ii,r*i) and (t 2 ,?’* 2 ) that satisfy the timelike separation 
condition |r *2 — r*i| < [^2 — iij, there exists a timelike geodesic of A42 that connects these two events. 

We have thus proven in this section that Af 2 is a causal domain. We recall that this has two important implications: 
first, that the world function a{x°‘,x^ ) is well-defined globally on Af 2 and second that the retarded Green function 
Q is likewise globally valid on Af 2 - 


IV. GEODESIC DISTANCE AND VAN VLECK DETERMINANT IN M 2 


In this section, we return to the wave equation on 4 — D space-time. We point out how the world function and 
Hadamard form of wave equations on M 2 have a role in the A — D analysis, and we calculate numerically two quantities 
of relevance. These quantities are the world function itself and the van Vleck determinant. For the calculation we 
apply some of the results obtained above for geodesics on A42, and we note that Theorem 1 provides the foundation 
on which these results are built: the calculations below are valid globally on Al 2 - 

By separating the angle variables in the retarded Green function in conformal Schwarzschild space-time in the usual 
way via a multipolar decomposition, we can write 

-| 00 

Gr{x,x') = —'^{2i + l)gi{x^,x^')Pi{cos'y), (4.1) 

£=0 


where Pg are Legendre polynomials and satisfy the PDE for a Green function on the 2-d conformal space: 


where 


{P-L^)Qi 


-A'k—52{ 



p = n2 + 


1 

4 



L = e + 


1 

2 ’ 


(4.2) 


(4.3) 
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and ^2 is the d’Alembertian operator of A42- In coordinates (t, r*) it takes the familiar form: 


I 

J.2 


02Q£ 


d'^Qt d’^Qt 


—47r(52(a:"^ — ), 


(4.4) 


where ^ is a Green function in Al 2 - 

The Hadamard form for the retarded Green function is an analytic expression which is only valid in a normal 
neighbourhood of a spacetime point. In (1 + l)-dimensions, the Hadamard form for the retarded Green function is 

(e.g., mmy- 

Q{x^,x^ ) = 6{t — t')9{—a)U{x^,x^ ). (4.5) 

The two-point function U{x^,x^ ) is a Riemann function for the wave equation |19j . and may be expanded in the 
so-called Hadamard series 


U{x^,x^') = Y,Uk{x^,x^')ay 

k—Q 


(4.6) 


Synge’s world function a satisfies 

VactV^ct = 2cr, (4.7) 

with the initial conditions \\Taxi^xCf{x,x') = 0 and b(x{x,x') = gAsix). The Hadamard coefficients 

Uk,k > 0, in the (1-1-1 (-dimensional space-time satisfy certain recurrence relations in the form of transport equations 
along the unique geodesic from x"^ to umiHiiiD]. In particular, in (1 -|- 1) dimensions, C/q is the square root of 
the so-called van Vleck determinant A = A{x^,x^ ) and it satisfies: 

Uo = cr^V^Go = (1 - □ 20 -)Go, lim Uq = 1. (4.8) 


Regularity at ct = 0 fixes constants of integration. This approach applies for any wave operator of the form P = 
□2 +v^\7A + w, where and w are respectively a vector and a scalar field on A42. Note however, that Uq is defined 
independently of the wave operator (as of course is a). 

In the rest of this section we numerically calculate Synge’s world function a and the van Vleck determinant A(x, x') 
in the 2-d conformal space A42 along the unique timelike geodesic joining two points in A42- Typically, in applications, 
one is concerned with causally separated points. For null separations, tr = 0 and so the calculation of a for time-like 
separations is the issue of most relevance. 


A. Calculation of a directly from the geodesic equations 


For timelike geodesics, the world function is given by cr = — where Ar is the difference between final and initial 
values of the proper time t along the unique timelike geodesic joining two timelike separated points (ti,ri), {t2,r2)- 
From now on we constrain ourselves to a subcritical {E < Ec) timelike (e = —I) geodesic which starts off at an 
initial point (ti, ri), bounces off the potential barrier at the periapsis r = r_|_ and then reaches a final point (t 2 , r 2 ), in 
coordinates {t,r). We note that this is the least straightforward case, and incorporates all technicalities that would 
be encountered in the other cases. 

The calculation of cr((ti,xi), (< 2 ,^ 2 )) proceeds in two steps. First, using (3.411, we can evaluate the proper time 


taken by a particle with energy E to fall from r = ri to r = r_|_. We do this for a range of values of E. Then, we 
(numerically) determine the (unique) value of E that corresponds to this particle departing from r = ri at coordinate 
time t = ti and arriving at r = r+ at time t = t 2 - The relevant value of At can then be identihed from the Hrst step. 

For convenience, we consider the case that r 2 = ri. Then, in order to calculate the corresponding value of the total 
prop er tim e interval Ar, we calculate Eq.( 3.43[ ) (where the dependence on xq = ri is in the variable if) - see Eqs.(3.46) 
and (3.47)) and multiply the result by two. In Eigj^a) we plot the total proper time interval Ar as a function of the 
energy for the case ri = 6M (which corresponds to the innermost stable circular orbit in Schwarzschild space-time). 

Similarly for the coordinate time, from Eq.( |3.28[ ) we have At = t 2 — ti = 2Ai(E). This integral can also be solved 
in terms of various elliptic functions but we do not write in the result as it is not particularly illuminating. We plot 
the result in Eigj^b). Eq.( |3.28 ) gives At = At{E) (for a fixed r+), so in order to find E = E{At) we look for a zero 
of this equation. We plot the resulting values of E as a function of At in FigQa). We can then introduce the values 
of E in Eq.( |3.43| in order to finally obtain the proper time interval At as a function of At for a given r 2 = ri. We 

plot the corrresponding a = — as a function of At in Eigj^ (b). 
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FIG. 3: Total interval of (a) proper time At from Eq.(3.43l and (b) coordinate time At from Eq.(3.28l, as functions of energy 
in the case r 2 = ri = 6M. 



FIG. 4: 
and (b) 


Various quantities as functions of At for ri = 6M: (a) E = E{At) (obtained by zeroizing ‘At —2Ai(i?)’ in Eq.(3.28l), 
(obtained from Eq.(3.43l using the values of E obtained as in plot (a)). 


(Ar) 

cr = 


B. Calculation of a and A from transport equations 

In this subsection we shall calculate Synge’s world function a{x, x') and the van Vleck determinant A{x, x') in M .2 
by numerically solving transport equations that they obey. 

The transport equation Eq.(4.8) obeyed by the van Vleck determinant may be simply written as 


dr 




(4.9) 


where r is an afhne parameter along the geodesic joining x and x' and = V^^V^cr. In order to obtain a transport 
equation for the derivatives of the world function, we combine Eqs.(3.11) and (5.2) in [20] together with: 


_ dQ B , -^A „.C/^D T^C „.Dr^A 

u Q B;C — — J - f 1 y B — ^ bdU Q Cj 

dr 


where is the tangent to the geodesic and Q^b = o'^b — 5^b- The following equation then follows: 


(4.10) 


dQ B DT^C Dy^A ^ t/ aA yaC I /aA \ t?A C D 

— J - =u Q BD — u I CdQ B - {Q cQ B + Q b) — tR CBDU U . 

dr T 


In A42 j the nonzero Christoffel symbols are 

r — M „„ QM^ — 5Mr + 


r(r - 2M)' 


— —- 
^ tt — 


pi ^ r-3M ^ y 

^ tr ^ rti 


r(r — 2M) 


(4.11) 


(4.12) 
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the nonzero Riemann tensor components are 


_ (r-2M)(r-6M) 

ttr — A — trt') 


Tit _ 

rtr — 


r — 6M 
r2(r- 2M) 


= -R 


rrt 7 


(4.13) 


and the Ricci scalar is equal to i? = 


12M 


- 2 . 


We can solve Eqs.(4.9) and ( 4.11[ ) simultaneously by using the code in [2T] in order to obtain r = T(r) and 
along a specific timelike geodesic in M 2 - We plo t thes e results in Figj^ We checked that the plot of 
r = T(r) agrees with the calculation of the proper time using Eq.(3.43) (with the initial value xq = ri being replaced 
by a value of r along the geodesic), thus serving as a validation of both methods. 

In Figj^we plot the final values (i.e., upon return to r = ri) of cr = —(Ar)^/2 and of Uq = A^/^ calculated as 
above for a geodesic in A42 for specific values of At and of ri = r 2 - These final values are the relevant ones for the 
expression for the Green function in Schwarzschild spacetime that we use in m- 




FIG. 5: (a) Proper time, r = r(r), and (b) square root of the van Vleck determinant, = A^^^(r), along a subcritical 

timelike geodesic in Al 2 . The geodesic has energy E ~ 0.1892071/M, it starts at r = ri = 6M, it reaches the periapsis r = r+ 
and it ends at r = ri. The results in the dashed blue curve (bottom-half curve in (a) and left-half curve in (b)) and in the 
straight green curve (top-half curve in (a) and right-half curve in (b)) respectively correspond to the trajectories r : ri —>■ r+ 
and r : r+ —>■ ri and they are obtained using the method described in Sec jlV B The curve in dashed red (overlapping with the 
dashed blue curve) in (a) is obtained using Eq.(3.43l (with the initial value xo = ri being replaced by a value of r along the 
geodesic). 


^ 1/2 




FIG. 6: Final values (i.e., at the return at r = ri) of: (a) the world function a = —(Ar)^/2 and (b) the square root of the 
van Vleck determinant, A^^^, along the only (subcritical and timelike) geodesic in M 2 as functions of At for ri = r 2 = 6M. 
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V. CONCLUSIONS 


The main motivation for this paper is ultimately the study of the self-force problem in General Relativity (although 
our results may also be useful for other problems, such as in the calculation of the response of a quantum “particle 
detector”). Approaches to this problem that rely on the study of Green functions in black hole spacetimes have 
yielded significant geometrical insights. Our particular focus (see e.g. m) is to determine global information on 


Green functions in black hole spacetimes. In the case of the exterior Schwarzschild spacetime, as seen in Section IV 
this is aided by a thorough understanding of Green functions on the spacetime M 2 - Thus M 2 is a space-time of 
physical relevance: it appears naturally in the wave equation in Schwarzschild space-time when applying a technically- 
convenient rescaling of the scalar field. Exploiting the natural occurrence of M 2 will be greatly aided by the geometrical 
property we have established here - i.e. that M 2 is a causal domain. 

On the one hand, such a proof is of theoretical interest since being a causal domain is a rather rare property of 
space-times. For example, black hole space-times, where closed null geodesics are a ubiquitous feature |22], are not 
causal domains. Likewise, plane-wave spacetimes generically contain pairs of points connected by multiple geodesics, 
and so are not causal domains [23] . De Sitter space-time contains pairs of points not connected by any geodesic [2l] . 
Thus geodesic convexity fails, and de Sitter space-time is not a causal domain. Anti-de Sitter shares this feature, and 
furthermore, there exist pairs of points p,q for which H J~{q) is non-empty and non-compact thus anti-de 

Sitter space-time fails both criteria required for a causal domain. (However, it should be noted that this failure to be 
a causal domain does not prevent one from constructing global Green functions on de Sitter space-time |7|.). 

From a mathematical point of view, the property of being a causal domain is unstable - at least in the case of Al 2 - 
To see this, consider perturbing the metric so that the effective potential shown in Figure has a local minimum. 
This local minimum would act as a potential well, trapping particles in a particular range of energies. Two particles 
with different energies, both confined to this potential well, will meet repeatedly, providing examples of pairs of points 
connected by multiple geodesics. The perturbation of the potential, and so of the metric tensor (and curvature tensor) 
may be arbitrarily small, and so a small perturbation of the space-time violates its causal domain nature. 

On the other hand, the property of being a causal domain is also of practical interest. The reason is that, being a 
causal domain implies that the world function is defined globally and that the analytic Hadamard form for the retarded 
Green function is valid for any arbitrary pair of points in M 2 - This feature is exploited in a separate study m, 
where we derive global properties of the Green function in Schwarzschild space-time by expressing it in terms of the 
world function and the Hadamard coefficients (such as the van Vleck determinant) in M 2 - 

Finally, we note that the focus on M 2 allows us to understand more clearly the global structure of Schwarzschild 
space-time in the following way. Consider null geodesics on Schwarzschild space-time. By conformal invariance of 
null geodesics, these are also null geodesics on the conformal Schwarzschild spacetime. The motion then projects onto 
geodesic motion on S 2 and geodesic motion on M 2 - Suppose that there are multiple null geodesics connecting the 
4 — D points a:“ = and a:“ = (^ 2 ;''’ 2 : ^ 2 ; </> 2 )■ Then the causal domain property of M 2 (in particular 

its geodesic convexity) implies that the projection of these two geodesics onto M 2 yields the same geodesic in M 2 
(but parametrized by different affine parameters). The projections into S 2 yield geodesics on S 2 that can differ only 
in terms of the particular great-circular path that they follow. Thus null geodesics that connect the same points of 
Schwarzschild space-time trace the same path in t — r (or t — r^) space. As exploited in m, this enables us to define 
a globally valid 2-point function on conformal Schwarzschild space-time that provides a global version of the world 
function. 
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